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TECHNICAL NOTE 3303 


TURBULENT -HEAT -TRANSFER MEASU REMEN TS AT 
A MACH NUMBER OF 3*03 
By Maurice J. Brevoort and Bernard Rashis 


SUMMARY 


Turbulent -heat -transfer measurements were obtained through use of a 
three-dimensional axially symmetric nozzle which consists of an inner 
shaped plug and an outer cylindrical sleeve. Measurements were taken 
along the outer sleeve and gave flat -plate results that are free from 
wall interference and corner effects for a Mach number of 3*03 and for 

a Reynolds number range of 5*6 X 10^ to 6.5 X 10^. The heat -transfer- 
coefficient results are in good agreement with theoretical analyses and 
the recovery-factor results are in good agreement with extrapolations of 
lower Reynolds number data. 


INTRODUCE ION 


The design of supersonic aircraft and missiles requires engineering 
information about heat -transfer coefficients and recovery factors that 
extend over a wide range of Reynolds numbers . Each of the existing experi- 
mental techniques for studying heat -transfer characteristics appears limited 
in various ways. A test setup which utilizes a flat plate for the testing 
surface is subject to corner and edge effects and also shock interference 
if high Reynolds numbers are obtained by extending the plate. 

One existing technique which avoids these effects is that of the 
hollow cylinder mounted centrally in the test section. However, the 
apparatus is subject to the tube choking at low Mach numbers and deviation 
from flat-plate conditions (the boundary- layer thickness approaching a 
sizable fraction of the cylinder radius) when the Reynolds numbers are 
low. In addition, there are differences in the heat -transfer character- 
istics for the inside and outside surfaces of the cylinders. 

The purpose of this paper was to present results obtained by using 
a technique that would avoid the above-mentioned effects and permit sim- 
ple and accurate measurements of local heat-transfer coefficients and 
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recovery factors over an extended range of both Reynolds and Mach numbers . 
The development- of-this technique for Mach number 3-03 and the resultB of 
preliminary heat-transfer measurements for. the cooling of the surface by 
the airstream are the subject of this paper. 


SYMBOLS 


c 


g 

h 

k 


m 

Nu 

Pr 


R 


St 

T 


eq. 


T 

T 

T 


stag 

v 

CO 


specific heat of-the sleeve material, Btu/lb-°R 
specific heat of the air at constant pressure, Btu/lb-°R 


acceleration due to gravity, ft/sec^ 
heat -transfer coefficient, Btu/ft^-see-°R 
heat conductivity, Btu/ft-sec-°R 
mass of- the sleeve material, pounds/ft 2 
Nusselt - ' number, hx/k 
Prandtl number, pCp^/k 

Reynolds number, pVx/|i 


recovery factor, 


T - T 
eq. 

T , - T m 

stag “ 


Stanton' number," Ku/RPr 


effective stream air temperature at the wall, seme tempera- 
ture which gives a thermal potential which is independent 
of -conductance h, °R 

stagnation temperature, °R 

inside surface temperature' - of L the sleeve, °R 
free-stream temperature-,- -°R 
free-stream velocity, ft/sec 


V 



HACA TN 3303 


3 


x longitudinal distance along sleeve, ft 

p free-stream density of the air, slugs/ft^ 

p dynamic viscosity coefficient, Tb-sec/ft 2 

dT/dt time derivative of average vail temperature for a particular 

x-station, °R/sec 

APPARATUS AND METHOD 


The apparatus consisted of a three-dimensional axially symmetric 
nozzle which was directly connected to the settling chamber of a cold-air 
blow-down jet. The nozzle consists of a shaped plug and an outer sleeve 
of 8-inch extra-heavy seamless carbon-steel pipe which was machined inside 
and outside to $ wall thickness of O .388 inch. A detailed drawing of the 
apparatus is shown in figure 1. The ordinates of the plug are given in 
table I, The plug was designed by using three-dimensional characteristics 
(see ref . 1) . 

Details of the installation of the thermocouples and static-pressure 
orifices are shown in figures 2 and 5- The thermocouples are located 
0.060 inch from the inner surface of the sleeve. The wires are number- 
30-gage copper- const ant an. (O.OIO inch in diameter). For these wires, 
conduction along the wire length is negligible. As indicated in figure 2, 
the thermocouples are in intimate contact with the steel wall so that 
thermal resistance at this junction is negligible. 

Figure 1 shows the location of the static -pres sure orifices. The 
Mach number distribution was measured both lengthwise and around the 
sleeve. The results are shown in figure 4. 

The temperature-recording equipment consisted of three synchronized 
high-speed Brown recorders, having a total of 36 switches, two of which 
were connected to the sett ling- chamber thermocouples, the rest being con- 
nected to the sleeve thermocouples. The stagnation-temperature variation 
with time for a settling-chamber pressure of 220 pounds per square inch 
gage is shown in figure 5 • Figure 6 shows the wall-temperature variation 
with time for stations 12 and 18 for a settling-chamber pressure of 
220 pounds per square inch gage. Test runs were made for settling-chamber 
pressures of 92, 142, 220, and 290 pounds per square inch gage. Excluding 
the first 20 seconds, the pressures were maintained constant for each run. 
The recorders were calibrated immediately before and after each run and 
are accurate to ±1° F. 
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In figure 7 there are plotted for various times during the non the 
curves of wall temperature against^ longitudinal distance, along the cylin- 
der. These curves were used to determine the rate of change of- the longi- 

tudlnal conduction k^-±. along the cylinder. Test results were taken 

dx 2 

only for the length of the cylinder for which k — - was zero. 

dx 2 

REDUCTION OF DATA 


Short test runs in which wall-temperature — time histories are 
obtained give the essential data from which heat-transfer coefficients 
and recovery factors may be determined. The reduction of these data 
requires a method in which the stagnation temperature and, accordingly, 
the equilibrium temperature may vary in an arbitrary manner. 

The Stanton number is calculated from 


St = 


h 

pvc pg 


and the he at - 1 ran s f e r coefficient is calculated from 


( 1 ) 


h = me 


dT/dt 

T - T 
A v - L eq 


( 2 ) 


By definition, the recovery factor is 




eq 


- T 


T 

stag 


(3) 


The method of- reducing the data is simply to select a recovery fac- 
tor and from equation ( 3 ) obtain T e ^. Substitute this T e( ^ into equa- 
tion ( 2 ) and obtain values of h for different heat flow rates. These 
values of h are then substituted into equation (l) . The true values 
of T e( ^, T| r , and St are -obtained when the St is constant with time 

(different heat flow rates). Figure 8 shows the curves used in evaluating 
the Stanton number and the recovery factor for station 12 for a settling- 
chamber pressure of 220 pounds per square inch gage. 
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The values of specific heat and specific gravity of the sleeve 
material were taken from reference 2. Figure 9 gives the variation of 
the specific heat of the sleeve material with temperature. The values 
of viscosity and Prandtl number ( 0 . 71 ) for air were also taken from 
reference 2 and were based upon the inside surface temperature of the 
sleeve. The value of T w chosen was at 120 seconds after starting. 
(This arbitrary choice of temperature can be made because the viscosity 
coefficient is practically a constant factor of the heat -conduct ion 
coefficient . ) 


ACCURACY 


The basic data consist of temperature-time histories at each thermo- 
couple location. The recorders were calibrated for synchronization imme- 
diately before and after each run and are accurate to +1° F. The thermo- 
couple wires are number-30-gage (0. 010-inch diameter) copper-constantan 
wires . Conduction along the length for these wires is negligible . The 
wires are insulated except at the junction where intimate contact with 
the sleeve keeps thermal resistance to a minimum. 

The sleeve is enclosed in a vacuum jacket to avoid free convection 
and, for the temperature range considered, radiation effects range up to 
2 percent. The product of the mass per unit area and the specific heat 
of the sleeve material as used in equation (2) is accurate to *5 percent. 

The value of dT/dt used in equation (2) is that of the actual 
thermocouple reading, whereas it should be the value associated with the 
average temperature in the thickness of the sleeve at each point. During 
the test runs, the inner-surface temperatures are lower than the average 
temperatures by approximately 1° to I. 5 0 at the beginning of the runs 
(40 seconds after starting) and become the same at the end of the runs 
(200 seconds). This effect results in the inner-surface temperature-time 
derivative being less than the .average temperature-time derivative by 
approximately 1 percent . The computed slopes when integrated reproduce 
the original curves to within ±0.5° F. 

The analysis of reference 3 was used to determine the effect of the 
sleeve diameter. This effect was computed to be less than 0.01 percent. 


RE SUITS AND DISCUSSION 


Figure 7 shows the variation of wall temperature with longitudinal 
distance along the cylinder for a sett ling- chamber pressure of 220 pounds 
per square inch gage. Over the test range, the wall temperatures are con- 
stant with x. In the relation 
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T w - T e q 

m and c are constant—and dl/dt is constant-because T v is constant. 
Therefore, if there is to be a variation of h with Reynolds number or 
x, T e q must vary with x. The value of T e q obtained for the run with 

a chamber pressure of 220 pounds per square inch gage, evaluated at 
120 seconds after starting, actually decreases approximately 2°. 

Figure 10 shows the variation of— local Ttfusselt numbers with local 
Reynolds numbers. The x used in evaluating these numbers was based 
upon an x that was considered zero at~the minimum. A single line 
faired through all the data points has a slope of approximately eight- 
tenths but, if lines are faired through the points for the individual 
runs, the slopes are slightly higher. It is reasonable to assume, how- 
ever, that the higher stagnation pressures would produce earlier transi- 
tion and hence the zero value of x would move upstream with increasing 
stagnation pressure. Hence, if" it— is assumed that the individual runB 
should merge into a continuous line (rather than a series of steps), x 
may be adjusted for x = 0 locations which would make all the runs coin- 
cident-.- This adjustment has been done in figure 11 by using values of 
x equal to zero at - 2.00, 4.25, 7*25, and 10.00 inches ahead of the mini- 
mum for the runs with chamber pressures of 92, 142, 220, and 290 pounds 
per square Inch gage, respectively. This systematic adjustment with 
pressure of the x = 0 location does not - affect the value of a faired 
curve through all the data^(figs. 10 and 11). 

The Nusselt numbers were found to vary from 4,125 to 32,500 for the 
Reynolds number range of 5*6 X 10^ to 6. 5 X 107 (based on adjusted x = 0 
locations). For comparison, the curves based upon the analyses of- Van 
Driest (ref. 4) and Eber (ref. 5) are shown. In the Van Driest analysis, 
T-^/Too was considered to be 3*0 (the average value of the test results). 

For comparison with these references, the data were computed by using the 
free-stream temperature to determine the density and the velocity. The 
wall temperature was used to determine the viscosity and Prandtl number. 
The two experimental values from V-2 data (ref. 6) are also determined 
in this manner. 

The results of this paper are approximately 10 percent higher than 
the Van Driest - analys i s and approximately 27 percent higher then the 
Eber analysis. The V-2 experimental results coincide with the results 
of this paper. 

Figure 12 shows the variation of-local recovery factors with local 
Reynolds numbers. The variation is from 0.887 at 5.6 X 10^ to 0.868 at 
6.5 X 107. The results are compared with an extrapolated curve which 
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represents the fairing of the data of reference 7* The data of refer- 
ence 7 are T° r Reynolds numbers below 6.75 X 10&. The results of this 
report and those of reference 7 agree to within *0.5 percent. Also 
included for comparison are the curves for recovery factor equal to 


Pr 1 /^ and Pr 1 / 2 . 
Prandtl number. 


The wall temperature was used to determine the 


CONCLUDING REMARKS 


It can be concluded that the method of testing and reducing the 
data gives re suit s that are in good agreement with established theoretical 
and experimental work. The reduction of the heat-transfer data requires 
adjusting the x = 0 locations (x is the longitudinal distance along 
sleeve) so that they are consistent when the Reynolds numbers are varied 
either by changing the pressure or the "x value" location. The varia- 
tion of Nusselts number with Reynolds number before and after adjustment 
of the x = 0 locations are essentially the same. 

The recovery-factor results are in good agreement with the extra- 
polated results of lower Reynolds number data. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 2k, 195k. 
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NOZZLE ORDINATES 

• • • * 


x, in. 

Y , in. 


-x, in # . 

y , in. 

Jo • 

* t 2 

3.532 


0 

3-532 

3.51^ 


.2, 

• 3.527 

.4 

3.4636 


•4 . 

3.512 

.6 

3.3905 


.6* * 

3.482 

.8 

3-3047 


.8 

3.444 

1.0 

3.2156 


1.0 

3.387 

1.6 

2.9612 


J-.5 

. 3.170 

2.0 

2.8036 


2.0 . • 

2.893 

2.6 

2.5841 


•2.5 

2.620 

3.0 

2.4482 * 


3-0 

2.380 

5.6 

2.2596 


3-5 

2.208 

4.o 

2.1436 


4.0 

2.100 

4.6 

1.9836 


4.5 

2.020 

5-0 

1.8856 


• 4.7 

2.000 

5.6 

6.0 

1.7509 

1.6687 


10.25 

2.000 

§ 



6 .G 

1.5565 




T.o 

1.4888* 




8.0 

1.3431 




10.0 

1.1493 

* 



12.320 

1*0764* 

• , 


• 

15.0 

1.0484 




20.0 

.$934 




25.0 

•*9384 




29.0 

•.8944 
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Figure 1.- Test arrangement 




Each wire wrapped and covered 
Both wires form one lead from 
each thermocouole . 



Sleeve (O.368 tbII thick) 
'' (steel) 


Figure 2.- Thermocouple installation. 
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Figure 5 .- Stagnation temperature variation with time for a settling 
chamber pressure of 220 pounds per square inch gage. 
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figure 8.- Stanton number as a function of time and recovery factors for 
station 12 for a settling-chamber pressure of 220 pounds per square 
inch gage. 
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Nusselt number 





Figure 10.- Variation of local Nusselt number with local Reynolds number 

for x = 0 at station 0. 
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Figure 12.- Variation of local recovery factor with local Reynolds 
number for corrected . x = 0 locations. Viscosity determined 
for wall temperatures. 





